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D . Using the multi-dimensional wave function formalism, we investigate the quantum state of a 

^"*^ ' scalar field inside a true vacuum bubble nucleated through false vacuum decay in flat spacetime. 

We developed a formalism which allows us a mode-by-mode analysis. To demonstrate its advantage, 
we describe in detail the evolution of the quantum state during the tunneling process in terms of 
individual mode functions and interpret the result in the language of particle creation. The spectrum 
of the created particles is examined based on quantum field theory in the Milne universe. 
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' I. INTRODUCTION 

Q\ 

Various phenomena associated with phase transitions at the early stage of the universe have been a subject of 
O" 1 great interest in cosmology for two decades. Cosmological first order phase transitions are caused by the decay of 
5^ , a metastable (false) vacuum in field theory, which produces true vacuum bubbles in the sea of false vacuum. A 
particular interest was taken in the false vacuum decay during an inflationary stage of the universe Jl]-|j| . 

Recently, several authors considered a possible inflationary universe scenario which realize Qo ~ 0.1 open universe 
|^-^|, in contrast to the standard inflationary universe which predicts £Iq = 1 flat universe Q). In this scenario, the 
bubble nucleation process plays an important role. That is, one assumes a scalar field with a potential like in the 
new inflationary scenario but with a high potential barrier before the slow rolling inflationary phase. If the nucleation 
rate is small, nucleated bubbles do not collide, and a homogeneous and isotropic universe with negative curvature is 
realized in one bubble thanks to the 0(4)-symmetry of the bounce solution. The second inflation starts inside a bubble 
due to the flat part of the potential. Then after the inflation, the universe gradually becomes curvature-dominated 
as it expands and the present universe with £7q ~ 0.1 can be explained. The important problem is the quantum state 
of the scalar field after bubble nucleation because its quantum fluctuations may be considered as the origin of the 
cosmic structure. 

As is well-known, the false vacuum decay rate can be calculated by means of the Euclidean path integral dominated 
by an 0(4)-symmetric bounce solution and the classical motion of the nucleated bubble is described by the analytic 
continuation of the bounce solution But studies of the quantum state after the false vacuum decay have not 

been done so much. Several pioneering works on this problem were done by Rubakov p0| , Kandrup JTT[ | , Vachaspati 
and Vilenkin [ fl2|| . Recently, we developed a formalism to investigate the quantum state after false vacuum decay 
based on the WKB wave function in a multi-dimensional tunneling system [|l3H16| , which was originally developed by 
Banks, Bender and Wu [l7|j , and Gervais and Sakita [T^ ]. 

For a simplified model of spatially homogeneous false vacuum decay discussed by Rubakov JhJ , we showed that our 
formalism leads to the same results |R|. As for the false vacuum decay due to an 0(4)-symmetric vacuum bubble, 
we also applied our formalism to a simple toy model and evaluated the expectation value of the energy momentum 
tensor inside the vacuum bubble The result indicates that the field excitation occurs inside the vacuum bubble, 
i.e., particle creation occurs there and it resembles a thermal state. 



In spite of these investigations, we have not yet obtained a clear understanding of the quantum state after bubble 
nucleation. Neither the spectrum of created particles, nor the relation between the efficiency of the particle creation 
and the shape of the field potential has not been made clear. A difficulty originates from the fact that the mass of the 
quantum field is spatially inhomogeneous and time-dependent due to the non-trivial background field configuration. 

In this paper, we tackle the issue by looking at the evolution of mode functions of a scalar field in a systematic 
manner and by investigating the spectrum of created particles, to obtain a better understanding of the quantum 
state after bubble nucleation. Our strategy is as follows. The background field configuration is inhomogeneous 
and time-dependent in the conventional Minkowski space-time coordinates. However, it keeps the Lorentz group 
(0(3, 1)) symmetry, which originates from the 0(4)-symmetry of the bounce solution. This implies that the system 
is homogeneous on the hyperbolic time slicing inside the future light cone within one bubble. This time slicing of 
the Minkowski spacetime is called the Milne universe. Quantum field theory in the Milne universe has been well 
investigated as an example of the quantization of a field in curved spacetime 19-22 , and the relation of the vacuum 
state between the Milne universe and the Minkowski spacetime is well understood for a massive field, at least for two- 
dimensional spacetime. Using this fact, we can define a natural vacuum state in the Milne universe that corresponds 
to the usual Minkowski vacuum. Then we can gain an insight into the quantum state of a field inside an expanding 
bubble in terms of the particle spectrum observed by the comoving observers in the Milne universe at late times. 

Following this strategy, we first carefully investigate the evolution of the mode functions through the tunneling 
process in two-dimensional spacetime, based on the wave functional formalism recently developed by us [ p"3| . We find 
that the resulting mode functions can be expressed in a very transparent manner and the created particle spectrum 
can be evaluated mode by mode. Then we extend the result to the case of four dimensional spacetime. This technical 
advantage will be very useful when we consider realistic cases of bubble nucleation. 

The paper is organized as follows. In section 2, we review our formalism to investigate the quantum state after 
bubble nucleation. In section 3, we investigate the quantum state inside a vacuum bubble in two-dimensional spacetime 
on the basis of quantum field theory in the Milne universe, and derive the formula for the spectrum of created particles. 
Then the result is extended to the four-dimensional case in section 4. In section 5, we apply our result to a simple 
thin- wall model which is similar to the one discussed by Rubakov |pxj| . The results are found to be similar to those 
obtained by Rubakov in the case of a spatially homogeneous decay model. Section 6 is devoted to summary and 
discussions. We use the units h — 1 and c = 1 and a bar (~) to denote the complex conjugate. 



II. REVIEW OF FORMALISM 



First, let us start reviewing our formalism to investigate the quantum state after false vacuum decay p3[ . It is 
based on the WKB wave function for a multi-dimensional tunneling system. For generality, we consider the (n + 1)- 
dimensional Minkowski spacetime. The conventional Minkowski spacetime coordinates are denoted by Jt, x) with 
x = (xi, ■ ■ ■ , x n ). We consider a system of two interacting real scalar fields a and 4> with the Lagrangian Jig] , 



C = C„ + Cd 



(2.1) 



where 



Co 



Ca = - / d 



Maf + V(a) 



(2.2) 
(2.3) 



Here the a field has a potential V(a) as shown in Figure 1, and decays from the false vacuum at <r(x) = up to the 
true vacuum at <r(x) = ut through bubble nucleation. 

In the lowest WKB picture, the decay from the false vacuum can be described by a bounce solution The bounce 
solution satisfies the classical Euclidean field equation, and has the maximum symmetry to minimize the Euclidean 
action. In our case, it has the 0(n + l)-symmetry, and depends only on r 2 +x 2 , where r (= it) is the time coordinate 
in the Euclidean space. Thus we can write the bounce solution as <tq(t 2 +x 2 ). The subsequent motion of the bubble 
after nucleation is given by the analytic continuation of the bounce solution to the Lorentzian time, which is given by 

CT (-i 2 +X 2 ). 

For simplicity, we neglect the quantum fluctuations of the a field and consider only the quantum state of the </> field 
which will be affected through the coupling with a as it decays from the false vacuum. Note that one may regard cf> 



2 



as the fluctuation of a itself, if desired (g6[. The wave function that describes the quantum state in the Euclidean 
region is written as |]l3| , |l4| 



where 



Se[<to] = f dr' fd n i 

J — OO J 



2 V 8t> 



+ V(a ) 



iPe[t, </>(•)] = cxp 



\J J d"xd n y0(x)^(x,y;r)0(y) 



(2.4) 

(2.5) 
(2.6) 



In Eq.(2.4), e -5 '^ ' ] is the lowest WKB part and gives the classical picture of the tunneling described by the bounce 
solution (To. The factor M(t) describes the next WKB order corrections to the normalization factor and is independent 
of (f>. The part is the wave functional that describes the state of the field <j> on the classical background of Co- 
Thus all the information of the quantum state is contained in the function f2e(x, y; t). It is expressed in terms of the 
Euclidean mode functions <?k('?~, x) as 



r» B (x,y;r) := J d n k 
where g^{r, x) satisfies the Euclidean field equation 

d 2 



rai <9g k (r,x) _ x 



Or 



8t 2 



E^-- 2 (-o(r 2 +x^)) 



i=i 



ffk(r,x) = 0, 



with the boundary condition at t — * -co as 



<? k -» e a "= T Y k (x), 



(2.7) 



(2.8) 



(2.9) 



where u>k '■— y k 2 + 77i 2 ((7f), up = <tq(t — > — oo), and lk(x) is a spatial harmonic function. This boundary condition 
is fixed by requiring that the system is set initially in a ground state in the false vacuum. The inverse of (^(r, x), 
t/k( r i x ) _1 is defined by 



J d n X 9k (T,x) flk ,(T,x)- 1 =5(k-k'). 



(2.10) 



The quantum state after false vacuum decay is described by the analytic continuation of the above wave functional 
into the Lorentzian region (r — > it) where the classical motion is allowed |13|] : 



where 



%o] 



dt'£ a [a (-t ,2 + X 2 )}, 



i>L[t,<f>(-)] = exp 



d"xd"y0(x)fi(x,y;t)0(y) 



(2.11) 

(2.12) 
(2.13) 



and ri(x, y;t) is expressed in terms of the mode functions u k (i,x), which is the analytic continuation of .g k (r, x) to 
the Lorentzian time, as 



n(x,y;t) / d»k 



1 /'^i, auk ( t ' x )r^T 



9* 



(2.14) 



Just as in the case of Eq.( |2~4f ), the first part e tS ^ in Eq.(2.11) describes the classical motion of the expanding bubble, 
and the second part ipL[t, </>(■)] does the quantum state of the <j> field on the classical background of the expanding 
bubble ct (x 2 - t 2 ). 
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It has been shown that the quant um state of the <j> field after false vacuum decay can be understood in the Heisenberg 
picture by using the mode fu ncti on u k (t, x) that satisfies the field equation for <j) on the classical background <7o |l3| , |l6| . 
As the boundary condition (|2.9|) does not fix the normalization of g k (r, x) and the function f2(x, y; r) is independent 



of it, we may assume that the mode function u k (t, x) satisfies the Klein-Gordon normalization (if it is not the case, 
we can take a linear combination of them to satisfy the normalization, i.e., v p — J d"kA(p, k)u k ), 

-i J d n x(u k (t, x)u k , (t, x) - ii k (t, x)u k , (t, x)) = <5(k-k'), (2.15) 

where a dot denotes the ^-differentiation. Introducing the creation and annihilation operators a k and a k , respectively, 
associated with these mode functions, the Heisenberg field operator <f>jj may be expanded as 

^ H = J d n k (a k u k {t, x) + 4u k (t,x)) . (2.16) 

Then the quantum state described by ipi,[t, (/>(■)] is found to be a "vacuum" state, defined by 

a k |^)=0 for any k, (2.17) 

in the Heisenberg picture. Note that this "vacuum" state is non-trivial in general. If there exists a natural choice 
of vacuum, say |0), the state is a squeezed state over |0), described by a non-trivial Bogoliubov transformation. 
Then it is possible to discuss the particle creation by evaluating the Bogoliubov coefficients between the two states. 



III. PARTICLE SPECTRUM 



In this section we first briefly review quantum field theory in the two-dimensional Milne universe |19|-|22j . Then we 
present a method to investigate the quantum state of a field inside a true vacuum bubble. Here we focus on the case 
of two-dimensional spacetime. The extension to the four-dimensional case will be given in section 4. 

The Milne universe is the hyperbolic time slicing of the Minkowski spacetime in the future light cone (see Fig. 2). 
Introducing the coordinates T and x m the two-dimensional Milne universe as 



t = Tcoshx, x = Tsinhx, (0 < T < oo, -co < \ < oo) 
the line element is written as 

ds 2 = -dT 2 + T 2 d X 2 - 
The field equation of a massive scalar field with mass M is 



d 2 



1 d 



i d 2 



dT 2 TdT T 2 d x 2 



4>(T,x) = 0- 



(3.1) 
(3.2) 

(3.3) 
(3.4) 
(3.5) 



where and are constants, and H„(z) and hL 2 \z) are the Hankel functions of the first and second kinds, 
respectively. Following the prescription of the second quantization of a field, we obtain 



Taking the mode expansion, 



<KT, X ) 



dp (p p (T)- 



2tt 



we easily find the solution as 



dp 



<Pp{T)- 



-=bp + tp p (T)—== 

' Ztt V27T 



hi 



(3.6) 
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where b p and bj, are the annihilation and creation operators, respectively, and \c p 2 ^ \ 2 — \c^ \ 2 = 1. The following choice 
of the coefficients have a special meaning. 



4^=0 cf = l. 



(3.7) 



Hereafter we regard b p and bp as those defined for this choice of the coefficients. The calculation of the propagator 
for the vacuum annihilated by b p leads to the same expression as that for the usual Minkowski vacuum |l9|,^o| , 



(01^(^)^(^)10) = ^H^(-M 2 (x^ - x'n 2 - ie). 



(3.8) 



(2) 

Thus H\ (MT) is the positive frequency mode function for the Minkowski vacuum in the Milne universe. 

For later convenience, we show this fact by directly relating the mode functions in the Milne universe and the usual 
Minkowski mode functions. By using the integral representation for the Hankel function 63], we find 



-it] cosh /C+iC sinh K. 



ip/2 



v - C 



(3.9) 



which holds under the condition, 

Im(?7 ± C) < 0. 

Setting ( — Mx, rj — Mt and JC — arcsinh(fc/M) in the above expression yields 

ip/2 



-iujfrt-\-ikx 



-iKp 



t + X 
t — X 



e ^ H < §(M(t 2 -x 2 ) 1 ' 2 ), 



(3.10) 



(3.11) 



where iv k = M cosh/C = Vk 2 + M 2 . If we use the coordinates (T, x) introduced in Eq. (|3.l|) , we have 



-iukt-\-ikx 



(2) 

The above formula means that H ip (MT) is indeed the positive frequency mode function for the Minkowski vacuum. 
Note that, because of the validity condition ( 3.10| ), the presence of a small negative imaginary part in t is understood, 
which corresponds to the familiar prescription. 

The negative frequency functions are also obtained by taking the complex conjugate of Eq.( 3.12] ), 



1 

2i 



dpe- lKp e ipx e np/2 H < >l ) (MT). 



(3.12) 



AuJkt—ikx 



l 

2i. 



dpe 



iKp -ipx„-Tp/2 rr(l) 



COT, 



where we have used the relation, 



e^/ 2 H^ 2) (MT) = e~^ 2 H^\MT). 



(3.13) 



(3.14) 



In contrast to the positive frequency mode function, Eq.(3.13) is valid on the upper half complex i-plane. 

Now we consider the mode functions describing the quantum state after bubble nucleation. Following the general 
discussion given in the previous section, we first need to solve the Euclidean field equation to find the Euclidean mode 
function 3fc(r, x). Introducing the 0(2)-symmetric coordinates £g and 6 in the Euclidean region as 



x = £ / ES,in8, r = — £gcos6>, (- 
the equation for the Euclidean mode function is written as 



19 1 92 V u 2 \\ 



d^E 2 ZEdtE tE 2 d9 2 



9k(T,x) = 0. 



(3.15) 



(3.16) 



The boundary condition is 3fe(r, x) 



a u> k T-ikx 



at r — ► — oo. 
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To find the solution, we first consider a solution of Eq.(3.16) in the form _F p (^)e pe . Then F p (£e) must obey 



d 2 



1 d 



d^E 2 £e di E £ 



fa -m 2 (M& 2 )) 



F p (Ce) = o. 



(3.17) 



In the limit ^ -> oo, we have Oq — ► up. Let M 2 = m 2 (<jE), then the asymptotic solution of -F p (£_e) which approaches 
zero at £e — ► oo is 



Fp&s) = K ip (M£ E ), 



(3.18) 



where K J/ (z) is the Modified Bessel function. On the other hand, from Eq.(3.13) with t — > —ir (r < 0), we find 



^uJkT — ikx 



dpe^ p e- pU K ip (MZ E ). 

j 

Thus the mode function g k (r, a;) can be expressed in the form, 

1 f°° 

g k (r,x) = - dpe^e-^F^), 
n J -00 



(3.19) 



(3.20) 



here F p (£e) is the solution to Eq.(3.17) with the condition F p (£e) — * Ki P {M^E) at £e — > oo. 

Once we obtain the Euclidean mode function, we continue it to the Lorentzian region by r — > it. This procedure 
can be readily performed in the region outside the light cone {x 2 — t 2 > 0) , 



u k (t,x) = - dpe iKp 
n J-oo 



t + x 
t — x 



-ip/2 



mi 



(3.21) 



where £ = V ' x 2 — t 2 . In order to obtain the mode function inside the future light cone, a bit more consideration is 
necessary. To keep the regularity of the mode function, we must perform the analy tic co ntinuation at an infinitesimally 
small Euclidean time before r = 0, r — > it — e(e > 0), or we must understand Eq.(3.21) with t replaced by t + ie [ pT[ . 
In terms of £, this implies £ = e~ ln / 2 T inside the future light cone. Then the mode function there is obtained as 
follows. In a sufficiently small neighborhood of £ = 0, we have m 2 (£) ~ M 2 = const.. Hence F p (£) can be expressed 
as 



F p (0 = A p K ip (M0 + B p (I ip (M0 + I- ip (M0) /2 
= A p K ip (MZ) + B p I ip (M£), 



(3.22) 



where A p — A p H — sinh7rpi?„. Note that, from the reality of FJ^e)i the coefficients A p and B p are real. Following 
the above prescription, the Modified Bessel functions in the above equation become 

K ip {MO -> K ip {Me-™'\t 2 - x 2 ) 1 ' 2 ) = ^ e -^ 2 H^(MT), 
I ip (MC) -> h P {Me~™/ 2 {t 2 - x 2 ) 1 ' 2 ) = l -e^' 2 \H\l\MT) + H {2) (MT)). 
Therefore we have the following expression for the mode function inside the future light cone at T ~ 0, 

Cp-e-^H^iMT) +T p e^' 2 H\l\MT)\ , 



u k {t,x) 



ilC P -ipx 



dpe l ^ p e 



where 



c p :- -A p 



cosh irp 



(3.23) 
(3.24) 

(3.25) 



(3.26) 



Taking the complex conjugate of Eq.( ft.2E ), we obtain 



G 



u k (t,x) 



dpe 



-iKp e ipx 



c p e^ 2 H%\MT) + d p e-^' 2 H^{MT) 



(3.27) 



where we have used Eq. (3.14). 

If we take the limi t of no er-dependence of mass, i.e., m 2 : (co(£e)) = M 2 — constant for all £e, we have c p = ir J2i 
and dp = 0, and Eq.( 3.27 ) reduces to 



u k (t,x) 



2i 



d V e- lKv e lvx e* vj2 H { ^ (MT) 



-iujk t-\-ikx 



(3.28) 



Thus the particle creation does not occur in this limit, as expected. 

In general, when the width of bubble wall cannot be neglected, the evolution of the classical background field inside 
the bubble becomes important. This effect leads to time variation of the mass of the <j> field, and additional particle 
creation may occur. In such a case, Eq.(3.27) merely gives the (asymptotic) initial condition of u k {t, x) at T — > 0. 



We now focus on the region inside the true vacuum bubble. Expanding the mode function as 



u k (t,x) = - 

7T 



dpe 



-iKp ipx 



(3.29) 



the field equation reduces to 



dr 2 + Tdf + T 2+m (ao( )} 



<P P (T) = 0. 



(3.30) 



thanks to the symmetry of the system. The initial condition for this equation is specified by Eq. ( 3.27| ) . We assume that 
m 2 (<Jo(T 2 )) — * fi 2 — constant at T — + oo, that is, ao(T 2 ) settles down to a finite value (presumably that corresponds 
to the true vacuum value oy). Then the asymptotic solution at T — > oo takes the form, 



/ \ 1 
u k (t,x) = - 



dpe' 



-ilCp e ipx 



(3.31) 



where c p and d p depend on the evolutionary history of m 2 (ao(T 2 )^ , 

In order to interpret the quantum state described by the mode functions obtained in Eq.(3.31) in the Heisenberg 
picture, we need to orthonormalize them. This is accomplished by taking the following linear combination of Uk(t, x): 



/oo 
dk\(q,k)u k {t,x), 
-OO 



(3.32) 



where 



X(q,k) = 



a iKq 



M cosh/C 



V 2 (i< 



Kl 2 ) 



(3.33) 



Then we have 



2tt 



2tt- 



(3.34) 



where 



\dq\ 



Pq--~- 



Cq\ 2 ~ \dq\ 



(3.35) 



<P q (T) 



Cq\ 2 



Kl 2 



--<Pq(T). 



(3.36) 
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One can easily check these mode functions satisfy the Klein-Gordon normalization on the T = constant hypersurface 
in the Milne universe 0] . 

The coefficients a q and (3 q are the Bogoliubov coefficients with respect to the natural vacuum in the Milne universe 
which corresponds to the usual Minkowski vacuum. The spectrum of created particles with respect to this natural 
vacuum is 



A, 



1 



_ 2 
Cq/dg - 1 



(3.37) 



When the particle creation after the bubble nucleation can be neglected, as is the case in the thin wall approximation, 
we have c„ = c„ 



^q and d q = d q . Then the spectrum of the particle can be expressed as 



f 



B, 



sinh 



B„ 



ttA„ 



(3.38) 



As clear from Eq.( 3.36| ), once we have the expression for ip p {T) in the bubble, it is trivially easy to obtain the 
orthonormalized mode functions v p (T, x)- In other words, the flat Minkowski coordinates play only an auxiliary role 
and we may focus on the evolution of a single mode function whose mode indices are specified with respect to (£e, &)- 
or (T, x)-coordinates. 



Let us summarize our result. First we solve Eq.(3.17) to find F p (£e) with the boundary condition that it decreases 
exponentially as £g — > oo. Then we analytically continue it to the interior of the future light cone F p (^ E ) — > <p p (T) by 
£e — > e~ I7r / 2 T and take its complex conjugate. The resulting function tp p (T) at T — > gives the initial condition of 
the mode function which describes the quantum state inside the bubble. If the evolution of the classical background 
field ao(T 2 ) can be neglected so that the mass of <f> is constant (= M.) inside the bubble, as in the case of a thin- 
wall bubble, <p p {T) at T — ► is the only information we need. By decomposing it into a linear combination of 

the asymptotic forms of e^ 2 H^(MT) and e^^H^ (MT) at T — > 0, the coefficients in front of them give the 
(unnormalized) Bogoliubov coefficients with respect to the natural vacuum of the Milne universe. The normalized 
Bogoliubov coefficients are then readily obtained as given in Eq.(3.35). If the evolution of the classical background 
cannot be neglected, as in the case of a thick-wall bubble, we must further solve the equation for the mode function 
with the initial condition specified by f P (T) at T — ► 0. Then at T — > oo where the mass of <fi settles down to a constant 

value (= fi), we again decompose <p p (T) into the two independent solutions e* p / 2 H^ (fJ,T) and e~ wp / 2 H^ (fiT) and 
read off the Bogoliubov coefficients with respect to the natural vacuum. We may then interpret the final quantum 
state in the language of particle creation due to bubble nucleation by considering the particle spectrum described by 
the Bogoliubov coefficients. 



IV. EXTENSION TO FOUR DIMENSION 



The analysis given in the last section can be extended to a four-dimensional system. We introduce the coordinates 
(T, x) in the four-dimensional Milne universe as 



Tsinhx, i = Tcoshx, (0 < T < oo, < x < °°) 



(4.1) 



where r and t are the radial and time coordinates in the Minkowski spacetime, respectively. The line element of the 
four-dimensional Milne universe is written as 



ds 2 = -dT 2 + T 2 (dx 2 + sinh 2 X^( 2 )): 
where dfl 2 2 ^ is the line element of the unit two-sphere. The field equation for a massive scalar field is 

d 2 



dT 2 ' TdT T 2 



0(T, x ,r») = o, 

where L 2 is the Laplacian operator on the unit three-dimensional spatial hyperboloid, 

L 2 



f 



9 ■ ( sinh 2 x S- ) + -V L ^> 



f 



sinh 2 % ®X \ dxl sinh 2 x 



(4.2) 



(4.3) 



(4.4) 
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and Lq is the Laplacian on the unit two-sphere. The eigenfunctions Y p i m of L 2 satisfy the equation 



and is in the form, 



where 



fplix) 



L 2 Y plm ( x ,n) = -(i + p 2 )Y plm ( x ,n), 

Y plm ( X ,n) = f pl ( X )Y lm (tl), 
T(ip + l + l) ] 



P i P -l/2( COsh X), 



T(ip) Vsmh^ vp - l l 



(4.5) 
(4.6) 

(4.7) 



^ 1/2 

and Yi m (Vt) is the spherical harmonics on the unit sphere, T(z) is the Gamma function and P ip _ 1 / 2 ( z ) * s the associated 
Legendre function of the first kind. The eigenfunctions Y p i m with < p < oo form a complete orthonormal set for 
square-integrable functions on the unit hyperboloid. Expanding the field operator in terms of Y p i m , we have 



fa(T,x,n)= / dp ]T 



,7rp/2 ip_ 



H%\MT) 



T 



plm 

+ h.c. 



(4.8) 



As shown in Appendix, the mode function e Kp l 2 (MT) /T turns out to be the positive frequency function with 
respect to the usual Minkowski vacuum, just the same as in the case of two dimension. 

Then we can repeat the procedure given in the previous section to find the mode function after bubble nucleation. 
Introducing the 0(4)-symmetric coordinates (£e,0) in the Euclidean region, 



r = £esin6', t=— £ecos( 



(0 < 6 < -, < £ E < oo) 



the equation for the Euclidean mode functions is written as 

d 2 3 d 1 



where 



S,e d£ E Ce 



1 d 



g k (r,x) = 0, 



J E — ~~T~2~ 

sin ( 



sin 



86 



sin 



(4.9) 



(4.10) 



(4.11) 



The boundary condition is <7k(r,x) — > e" fc7 ' ji(kr)Yi m (Q,) at r — > — oo. 
Using the formula ( |A9| ) in Appendix, g^(r,x) is expressed as 



5k (r,x) = (-1)' 



x sm 



p dp 



P,?,, 1 /^ (cosh AC) 
V sinh AC 
'sinhpfl F p (£ E ) 



d cos / sin M£e 



(4.12) 



where F p (^ E ) satisfies the same equation and the boundary condition as in the two-dimensional case, Eqs.(3.17) and 
(3.18). This fact leads to the following expression for the mode function inside the future light cone, 



i+i 



dp- 



T{ip + l + l)\ 2 P-^icoshlC) P7p$ (cosh*) 



MT 



\T{ip)\ 

0/2 rr(2) 



V sinh JC 



Vsinh x 



c p e^ 2 H^'(MT) + d p e-*P/ 2 H£ } (MT) 



(4.13) 



where c p and d p are the same as those defined in Eq.( |3.26 ). Thus the argument goes completely parallel to the 
two-dimensional case. 
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If the background field evolves inside the bubble, we must further solve Eq.(B.30) for the mode function with the 
initial condition at T — > given by Eq.( 4.13 ). The result will take the form, 



M k (i,x) = (-£) 



T( IP + l + l)\ 2 (coshJC) P+ l fi (coshx) 



dp 

o |r(ip)| 
1 



V sinh /C 



\/sinh x 



where we have assumed m 2 ((7o(T 2 )) — > ^i 2 = constant at T — > oo. 



To obtain the orthonormalized mode functions, instead of Eq.(3.33), we set X(q, k) as 

■ ^ (sinh/c) 3 / 2 r(?; g + i + i) p _,_ 1/2/ 



X(q,k) 



2 '{\c q \ 2 - Kl 2 ) cosh/c 



r(i?) 



ig-l/2 



(cosh/C), 



and take the linear combination, 



v q {T,x,ty= dfcA(g, fc)u k (i,x). 
Jo 



Using the orthonormality condition of the harmonic functions, 



/•OO 

/ d X sinh xPrp-i/i 2 ( cosh X^y-i/a ( c «sh x) 

«/ 



|r(*p)| 



|r(tj»+i + i) 



r 6(p-p')> 



(4.14) 



(4.15) 



(4.16) 



(4.17) 



we then find the mode functions which satisfy the Klein-Gordon normalization on the T — constant hypersurface in 
the Milne universe as 



T 



Pi 



T 



f q i(x)Yi m (n), 



(4.18) 



where the coefficients a q and f3 q are the same as t hose given by Eq.( 3.3li ). Thus the particle spectrum is also the 
same form as in the two-dimensional case, Eq.(3.37). 



V. SIMPLE MODEL 



In this section, applying the technique developed in the previous sections, we consider a simple model of the particle 
creation through the tunneling process of bubble nucleation. We assume a thin-wall bubble so that the mass of the 
cj) field is given by 



m 2 (cr (££ 2 )) = 

where R is the radius of the bubble when it nucleates. 



M 2 ; R < £, E < oo, 
< £e < R, 



m 2 ; 



The solution to Eq.(3.17) can be easily found by matching the solutions in both regions, 
where 



K ip (M£ E ); i?<£ B <oo, 
A p K ip (n£ E ) + BpIip(n^ E ); < £ E < R, 



A p = fiRIi P (jiR)Ki P (MR) - MRI ip (nR)K ip '(MR), 
B p = -fiRK ip '{(j,R)K ip (MR) + MRK ip (nR)K ip '(MR) 

and a prime denotes differentiation with respect to the argument. 

_ 2 



The particle spectrum is given by n p — 



Bp/nAp 



(5.1) 
(5.2) 



Roughly, the Modified Bessel functions asymptotically behave 



as Kip(z) oc e 7rp ^ 2 and Ii P (z) oc e np ^ 2 for p — > oo. Hence the particle spectrum always decreases exponentially for 
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large p as n p oc e~ 2lrp . This suppression of the spectrum at high momentum limit resembles that of a thermal state. 
The same feature has been found in a spatially homogeneous model of false vacuum decay, first discussed by Rubakov 
jlO| and subsequently by us jb|] . Note that the mode index p corresponds to a comoving wave number in the Milne 
universe and the scale p = 1 corresponds to the physical curvature scale on the hyperboloid. In what follows, we 
evaluate n p for several limiting cases. 

(1) M < Br 1 and fi < Rrh 

In this case, we can use the following asymptotic formulae for the Modified Bessel functions at z <C 1, 

( z /2yp 



lip(z) 



T(ip + 1)' 



R ip( z ) ~ o 777 



1 



(z/2) 



-lp 



(z/2) 



2 ismhirp\T(—ip + 1) T(ip+1 
After some algebra, we find the following expression for the particle spectrum, 

_ sin 2 (plog(/VM)) 

n P — . , 2 

Sinn np 



(5.3) 



(5.4) 



We see that n p decreases exponentially at large p, in accordance with the general discussion in the above. Due to 
the sinusoidal dependence, it has zeros in a period of ir/ \og{ji/M). In the limit p — > 0, n p approaches a finite value, 



log( M /M)/7r 



Figure 3 shows the spectrum for the case MR — 0.02 and fj,R = 0.01. 



(2) M > R- 1 and fj, < Br 1 : 

In this case, if we focus on the region p < 1, the asymptotic formulae for the Modified Bessel functions at z ^> 1 
can be used, 



hp(z) 



1 



/27TZ 



K lp {z) ~ J—e z , 



for z > 1. 



(5.5) 



Then n p is evaluated as 



1 



4 sinh irp 



T(ip) /fiR\- 2ip 
T(-ip) \~Y) 



(5.6) 



In the limit p — > 0, we find n p — > (7 + log(/ii?/2)) /n , where 7 = 0.577 • • • is the Euler constant. Figure 4 shows n p 
as a function of p for the case MR = 20 and fj,R = 0.01. 

(3) M > iT 1 and /1 > Rrh 

Using the asymptotic formulae ( |5.5| ), we easily find n p in the region p ^ 1, 



M + M 



(5.7) 



(4) M < i?- 1 and fi > ij" 1 : 



In this case, using Eqs.( |5.5| ) and (5.5), we obtain n p at p < 1 as 



(5.8) 



The first two cases for which fiR <C 1, the spectrum shows a similar feature, that is, it is constant of order unity at 
p < 1, and decreases exponentially at p > 1. On the other hand, in the latter two cases when ijlR 1, the particle 
creation is always suppressed exponentially by a factor of 0(e~ ifJlR ). This feature is the same as those obtained 
previously in the spatially homogeneous decay model of the false vacuum if we identify the Euclidean tunneling 
"time" in the homogeneous model with the radius of the bubble wall 1 1 l^Jl^,^]. 



11 



VI. DISCUSSION 



In this paper, we have investigated the quantum state of a scalar field inside a vacuum bubble nucleated in the 
process of false vacuum decay, focusing on the aspect of particle creation. We have used the wave functional formalism 
p3[ to investigate the quantum state. First, we have carefully analyzed the case of two-dimensional spacetime and 
presented a prescription to obtain the mode functions inside a future light cone of a bubble and an interpretation 
of the resulting quantum state in the particle creation picture. Because of the existence of a natural vacuum in 
the Milne universe, which corresponds to the conventional Minkowski vacuum, we have been able to introduce the 
concept of "particle" in a natural way. The spectrum of created particles follows immediately by reading off the 
Bogoliubov coefficients relative to the natural vacuum from the final form of the mode functions. We have shown 
that the Bogoliubov coefficients are diagonal in the mode indices so that a mode-by-mode analysis of the quantum 
state is possible. Then we have extended our method to four-dimensional spacetime and shown that it applies in 
essentially the same manner as in the case of two dimension. This technical advantage will be very useful when we 
consider realistic models, such as the case of a thick-wall bubble or the case of quantum fluctuations of the tunneling 



field itself 1 26 1 



We have applied the technique developed here to a simple thin- wall model and found the following. (1) The 
spectrum of created particles is always suppressed exponentially at large wavenumber with the characteristic scale of 
the spatial curvature on the hyperbolic time slicing inside the light cone (i.e., the cosmic time slicing in the Milne 
universe). (2) If the field becomes massive inside the true vacuum compared to the inverse of the wall radius, i.e., 
[i ^> R -1 , the particle creation is always suppressed by an exponential factor of 0(e~ 4fJ - R ). 

In this paper, we have assumed the spacetime to be flat and only considered a massive field inside the true vacuum 
bubble. However, it is of great interest to develop a formalism by relaxing these restrictions. In the case of a massive 
field, the quantum state corresponding to the Minkowski vacuum state can be written as a pure state in the Milne 
universe. But it will not be the case for a masslcss field ^7|. This is because the massless modes propagate always 
along null directions and a half of them escape to future null infinity without intersecting the boundary light cone 
of the Milne universe. Thus the massless Minkowski mode functions can never be expressed in terms of those in the 
Milne universe alone. Apparently a more careful analysis is necessary for the massless case. 

Extension to the case including gravity is most interesting, especially in connection with the open universe inflation 
^,[1 . As pointed out in the Introduction, if we consider a one-bubble inflationary universe model to realize a negative 
curvature universe, the spacetime will be a de Sitter universe and the quantum fluctuations of the inflation field inside 
a nucleated bubble will give rise to the cosmological density perturbations. Hence, if it is possible to include the 
effect of gravity in our approach, we will be able to determine the spectrum of the density perturbations in the open 
universe inflationary scenario. This is now under study ||. 
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APPENDIX: 

In this Appendix, we derive the formulae which relate the mode functions in the Milne universe and those of the 
usual Minkowski spacetime in four-dimensional spacetime. As we use the spherical coordinates in both cases, we omit 
the spherical harmonics from the mode functions. We start from the following formula which may be derived from 



Eq.(3.11) or (3.12) 



ke-^j (kr) = - r dp sin^^e^ 2 gg)(Mr) . (Al) 
Jo smh x T 

This gives the relation between the I = Milne and Minkowski mode functions. 

We extend the above formula to the case of I > 1. By mathematical induction we can show that 
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Using this fact, we operate (sinhx) (d/dcoshvj' on both sides of Eq.(Al) and integrates both sides of the equation 
with respect to u>k, we find 

where 

p(jj k pui i>(jJ2 rwi—\ 

a p k = I duii / du>2 / dui3 ■ ■ ■ / du>i smpK.(u>i) 
Jai Jai Jai Jai 

dw 2 (wfc-cj 2 ) / dw 3 ■ ■ ■ / dun sinplC(uii) 

AI JAI JAI 

u k ((),—() V -1 

**i „ A, sinp/C^). (A4) 

Noting that cj^ = Mcosh/C, if we change the integration variable u)\ to arccosh(wi/M), we see that this is just an 
integral representation of the associated Legendre function pa] , 



= M l ^-p (sinh/C)'+ 1 / 2 J P; i r i 1 / 2 2 (cosh/C). (A5) 



Hence we have 



roo r~ 

l ji(kr) = -i l pdp^j-- 



7T PJ-ifr(<xxh£) 



ip-l/ 



V sinh JC 



XSmh K^sh^J sinhT Mr ■ (A6) 

Using the formula, 

^. J' Vgg ■ til JE&ail. ilTS^f^, (A7) 

Vdcoshx/ sinhx P V 2 |r(ip)| V sm hx 
which can be derived by using the differential relation for the associated Legendre functions, we find 

|2 p-l-l/2. 



e-^Mkr) = / dp 



2 



2 Jo \T(ip)\ Vsinh/C 

Pr l ~i\C 2 (cosh y) , H (2] (MT) 

This result shows that e 7rp ^ 2 H^\MT)/T is indeed the positive frequency function in the Milne universe which 
describes the Minkowski vacuum, as in the tw o-dimensional case. 

By taking the complex conjugate of Eq.(A6), and continuing it to the Euclidean time t — > —it (t < 0), we obtain 

e^ji(kr) = (-l)\- pdp-- ip - 1/2[ ' 



7T Jo Vsinh K 

x sin < e ( -±- V ^ Be^eI . ( A9 ) 

V d cos J sin m^E 



This equation corresponds to Eq.(3.19) in the two-dimensional case and gives the expression for the Euclidean mode 
function in the four-dimensional case. 
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FIGURE CAPTION 
Fig.l. An illustration of a potential for the tunneling field a. 
Fig. 2. A sketch of the coordinates describing the Milne universe. 

Fig. 3. The number spectrum of created particles as a function of the wave number p for the case MR = 0.02 and 
fiR = 0.01. 

Fig.4. The same as Figure 3, but for MR = 20 and fiR = 0.01. 
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